On Gravity

by Jack Pickett - London & Cornwall - October / November 2025

Introducing a single, universal gravitational law...

geff=GMr2eκr g_{\text{eff}}=\frac{GM}{r^{2}}\mathrm{e}^{\kappa r}

Gravity is not only determined by mass and distance. It also depends on how matter is distributed in any given situation. The term κ measures how the local density environment influences the strength of gravity.

Matter bends the space around it — but how much it bends depends on how that matter is distributed.

Consider a kitten on a mattress. It will make no visible indentation on the mattress.
Now consider 1000 kittens all arranged in a grid on the mattress: still no visible indentations.

Now move more kittens into the center of the mattress and, gradually, an indentation will form.
Furthermore, if we swap the mattress for actual spacetime, and add a dense enough region of kittens, the curve becomes so deep that light can't escape and we are left with a black hole! (And kitten spaghetti...)

κ\kappa

In Newtonian gravity the potential follows a strict inverse–square form. This assumes a uniform mass distribution. Real systems contain gradients, shear, and density variations that alter the effective curvature.

A small modification to the curvature term in the action,

R    ReαRR \;\rightarrow\; R\,e^{\alpha R}

introduces an exponential correction to the weak–field potential. The resulting effective potential takes the form

Φeff(r)=GMreκr\Phi_{\rm eff}(r) = -\dfrac{GM}{r}\,e^{\kappa r}

and the corresponding acceleration becomes

geff(r)=GMr2eκrg_{\rm eff}(r) = \dfrac{GM}{r^{2}}\,e^{\kappa r}

The parameter κ enters as the weak–field imprint of this geometric modification. It encodes how local structure modifies the effective curvature.

κ(r)\kappa(r)

The geometric origin implies that κ depends on the local environment. Observationally, the dominant contributions arise from background curvature, velocity shear, and density.

κ=κ0  +  kv(v/r1012s1)3(ρρ0)1/2 \kappa = \kappa_{0} \;+\; k_{v}\, \left(\frac{\partial v / \partial r}{10^{-12}\,\mathrm{s}^{-1}}\right)^{3} \left(\frac{\rho}{\rho_{0}}\right)^{1/2}
κ₀ background curvature
kᵥ shear–response coefficient
∂v/∂r local velocity gradient
ρ density relative to ρ₀

κ increases in regions with strong shear or enhanced density, and decreases in smooth or diffuse environments. This produces the observed variation in gravitational behaviour across galaxies, clusters, and large–scale structure.

Geometric Algebra View (Optional)

For a deeper geometric intuition, consider κ in Clifford algebra Cl(1,3) — the language of spacetime rotations and oriented areas. In this picture, curvature generated by density gradients is represented by a bivector, an oriented plane element.

B=ereρ B = \mathbf{e}_r \wedge \mathbf{e}_\rho

where:

eᵣ = unit radial direction (along r)
eᵨ = unit density gradient, ∇ρ / |∇ρ|

This bivector B defines the plane in which radial paths react to structure. Its magnitude |B| measures how strongly matter clumping twists or redirects those paths. κ can be interpreted as an effective scalar built from |B|, encoding how local structure modifies the gravitational field.

This complements the f(R) derivation used above: exponential behaviour in modified gravity emerges naturally from geometric “wedges” in the Ricci curvature. Teleparallel analogues such as f(T) = T exp(βT) offer a torsion-based formulation where T ∼ |B|² links directly to the same bivector structure (Nojiri 2007; Farrugia 2016).

Vera Rubin stars

When astronomers calculated how fast stars should orbit in a galaxy, they used the standard intuition that stars near the center should orbit fast, and stars farther out should orbit much slower, because they are farther from most of the galaxy’s central mass. However Vera Rubin's observations contradicted this: the stars at the edges were not slowing down. They were moving just as fast as the stars near the center. In many galaxies, they move about three times faster than both Newton & Einstein predict.

Since κ adjusts gravity based on how matter is distributed, we can apply it directly to a real galaxy to see whether it reproduces the observed rotation speed:

Andromeda (M31) observed:  250 km s1 \textbf{Andromeda (M31) observed:}\approx\;250\ \text{km s}^{-1}
Newton predicts:vN=GMr \textbf{Newton predicts:}\qquad v_N=\sqrt{\frac{GM}{r}}
v_N ≈ sqrt(6.674e-11 * 2.0e41 / 8.0e20) ≈ 1.29e5 m/s ≈ 129 km/s 🛑
With curvature response (κ):vκ=vNeκr/2 \textbf{With curvature response (}\kappa\textbf{):}\qquad v_\kappa=v_N\,e^{\kappa r/2}
v_κ = v_N * e^(κr/2) with κ ≈ 1.65e-21 → v_κ ≈ 129 km/s * e^((1.65e-21 * 8.0e20)/2) ≈ 250 km/s ✅

To compare theory with rotation data, we derive κ directly from observations by taking the Newtonian speed from "baryonic" mass, compare to the observed speed, and solve for κ.

Galaxyradius (m)Mass (kg)κ (m⁻¹)Newton predicts (m/s)v_model (m/s)v_obs (m/s)
Milky Way3.086e201.2e412.0196e-21κ ≈ (2 / 3.086e20) * ln(2.2e+5 / 1.61096e+5) ≈ 2.0196e-21 m^-1🛑 161.1 km/sv_N ≈ sqrt(6.674e-11 * 1.2e41 / 3.086e20) ≈ 1.61096e+5 m/s✅ 220 km/sv_model ≈ 1.61096e+5 * exp((2.0196e-21 * 3.086e20)/2) ≈ 2.2e+5 m/s220 km/s
NGC 31989.26e201.0e411.43524e-21κ ≈ (2 / 9.26e20) * ln(1.65e+5 / 8.48961e+4) ≈ 1.43524e-21 m^-1🛑 84.9 km/sv_N ≈ sqrt(6.674e-11 * 1.0e41 / 9.26e20) ≈ 8.48961e+4 m/s✅ 165 km/sv_model ≈ 8.48961e+4 * exp((1.43524e-21 * 9.26e20)/2) ≈ 1.65e+5 m/s165 km/s
NGC 24033.086e202.0e404.66074e-21κ ≈ (2 / 3.086e20) * ln(1.35e+5 / 6.57673e+4) ≈ 4.66074e-21 m^-1🛑 65.77 km/sv_N ≈ sqrt(6.674e-11 * 2.0e40 / 3.086e20) ≈ 6.57673e+4 m/s✅ 135 km/sv_model ≈ 6.57673e+4 * exp((4.66074e-21 * 3.086e20)/2) ≈ 1.35e+5 m/s135 km/s
NGC 29033.703e206.0e403.53239e-21κ ≈ (2 / 3.703e20) * ln(2e+5 / 1.0399e+5) ≈ 3.53239e-21 m^-1🛑 104 km/sv_N ≈ sqrt(6.674e-11 * 6.0e40 / 3.703e20) ≈ 1.0399e+5 m/s✅ 200 km/sv_model ≈ 1.0399e+5 * exp((3.53239e-21 * 3.703e20)/2) ≈ 2e+5 m/s200 km/s
NGC 9254.63e206.0e409.17251e-22κ ≈ (2 / 4.63e20) * ln(1.15e+5 / 9.2999e+4) ≈ 9.17251e-22 m^-1🛑 93 km/sv_N ≈ sqrt(6.674e-11 * 6.0e40 / 4.63e20) ≈ 9.2999e+4 m/s✅ 115 km/sv_model ≈ 9.2999e+4 * exp((9.17251e-22 * 4.63e20)/2) ≈ 1.15e+5 m/s115 km/s
NGC 5055 (M63)8.95e203.0e415.92716e-22κ ≈ (2 / 8.95e20) * ln(1.95e+5 / 1.49569e+5) ≈ 5.92716e-22 m^-1🛑 149.6 km/sv_N ≈ sqrt(6.674e-11 * 3.0e41 / 8.95e20) ≈ 1.49569e+5 m/s✅ 195 km/sv_model ≈ 1.49569e+5 * exp((5.92716e-22 * 8.95e20)/2) ≈ 1.95e+5 m/s195 km/s
NGC 73311.08e214.0e417.83305e-22κ ≈ (2 / 1.08e21) * ln(2.4e+5 / 1.57221e+5) ≈ 7.83305e-22 m^-1🛑 157.2 km/sv_N ≈ sqrt(6.674e-11 * 4.0e41 / 1.08e21) ≈ 1.57221e+5 m/s✅ 240 km/sv_model ≈ 1.57221e+5 * exp((7.83305e-22 * 1.08e21)/2) ≈ 2.4e+5 m/s240 km/s
NGC 69464.32e201.3e418.42427e-22κ ≈ (2 / 4.32e20) * ln(1.7e+5 / 1.41717e+5) ≈ 8.42427e-22 m^-1🛑 141.7 km/sv_N ≈ sqrt(6.674e-11 * 1.3e41 / 4.32e20) ≈ 1.41717e+5 m/s✅ 170 km/sv_model ≈ 1.41717e+5 * exp((8.42427e-22 * 4.32e20)/2) ≈ 1.7e+5 m/s170 km/s
NGC 77931.85e208.0e396.16275e-21κ ≈ (2 / 1.85e20) * ln(9.5e+4 / 5.3722e+4) ≈ 6.16275e-21 m^-1🛑 53.72 km/sv_N ≈ sqrt(6.674e-11 * 8.0e39 / 1.85e20) ≈ 5.3722e+4 m/s✅ 95 km/sv_model ≈ 5.3722e+4 * exp((6.16275e-21 * 1.85e20)/2) ≈ 9.5e+4 m/s95 km/s
IC 25742.16e203.0e397.0226e-21κ ≈ (2 / 2.16e20) * ln(6.5e+4 / 3.04458e+4) ≈ 7.0226e-21 m^-1🛑 30.45 km/sv_N ≈ sqrt(6.674e-11 * 3.0e39 / 2.16e20) ≈ 3.04458e+4 m/s✅ 65 km/sv_model ≈ 3.04458e+4 * exp((7.0226e-21 * 2.16e20)/2) ≈ 6.5e+4 m/s65 km/s
DDO 1541.85e201.0e391.0464e-20κ ≈ (2 / 1.85e20) * ln(5e+4 / 1.89936e+4) ≈ 1.0464e-20 m^-1🛑 18.99 km/sv_N ≈ sqrt(6.674e-11 * 1.0e39 / 1.85e20) ≈ 1.89936e+4 m/s✅ 50 km/sv_model ≈ 1.89936e+4 * exp((1.0464e-20 * 1.85e20)/2) ≈ 5e+4 m/s50 km/s

TLDR: considering density distribution seems to matter. (dark matter...)

Gravitational Lensing

The next question is whether this same curvature term applies to light as well as mass. Gravitational lensing allows us to test that directly by comparing the bending of light predicted from observed mass to the bending we actually observe.

In galaxy rotation, orbital velocity depends on the square root of the gravitational potential. This means the κ effect shows up as a factor of exp(κ·r / 2). In gravitational lensing, the bending of light depends on the potential directly, not its square root. So the same κ shows up as exp(κ·b / 2), where b is the light’s closest approach to the mass.

αeff(b)=(4GMc2b)eκb/2 \alpha_{\text{eff}}(b) = \left(\frac{4 G M}{c^{2} b}\right) \mathrm{e}^{\kappa b / 2}
Same k - different observables
LensM (kg)b (m)α_GR (arcsec)κ (m⁻¹)e^(κ b/2)α_model (arcsec)α_obs (arcsec)
Abell 1689 (cluster)2.0e453.0e21🛑 408.45″
α_GR = 4GM/(c²b) → 0.001980220393 rad
-1.47047e-211.10173e-1✅ 45″
α_model = α_GR · e^(κb/2) → 0.000218166156 rad
45″
Bullet Cluster 1E 0657-5582.0e454.5e21🛑 272.3″
α_GR = 4GM/(c²b) → 0.001320146929 rad
-1.01098e-211.02828e-1✅ 28″
α_model = α_GR · e^(κb/2) → 0.000135747831 rad
28″
MACS J1149.5+2223 (cluster)1.0e453.6e21🛑 170.187″
α_GR = 4GM/(c²b) → 0.000825091830 rad
-1.13659e-211.29269e-1✅ 22″
α_model = α_GR · e^(κb/2) → 0.000106659010 rad
22″
SDSS J1004+4112 (quad QSO, cluster-scale lens)3.0e446.5e20🛑 282.773″
α_GR = 4GM/(c²b) → 0.001370921811 rad
-9.24796e-214.95097e-2✅ 14″
α_model = α_GR · e^(κb/2) → 0.000067873915 rad
14″

Collisions

During high-velocity cluster collisions, gas clouds experience shock compression and strong velocity shear, raising κ temporarily:

κ=κbase+κcoll \kappa = \kappa_{\text{base}}+\kappa_{\text{coll}}

where

κcoll=kv ⁣(vrel1012 s1) ⁣3(ρρ0) ⁣1/2 \kappa_{\text{coll}} = k_v\!\left(\frac{\nabla v_{\text{rel}}}{10^{-12}\ \mathrm{s}^{-1}}\right)^{\!3} \left(\frac{\rho}{\rho_0}\right)^{\!1/2}
kv5×1026 m1,  ρ0=1600 kgm3 \quad k_v \approx 5\times 10^{-26}\ \mathrm{m}^{-1},\ \ \rho_0=1600\ \mathrm{kg\,m^{-3}}

Gravitational lensing depends on the gravitational potential and increased κ multiplies the bending angle. As the shock and shear dissipate, κ_coll → 0 and the lensing map recenters naturally.

The lensing region shifts — appearing heavier — but "extra mass" is not needed when described as extra weight.

Bullet Cluster — Collision Shift

cluster separation: 160.0 px
κ_base = 7e⁻²¹ m⁻¹
κ_coll(t) = 1.58e-5
κ_total = κ_base + κ_coll
lensing amplification α_model / α_GR ≈ 1.00
apparent lensing center shift: 0.0 px

As the clusters pass through each other, the regions of strongest curvature shift — not because new mass appears, but because the collision briefly increases the weight of space itself.

Local Group — Basin Map

This map shows the gravitational potential of the Local Group as a continuous basin, where the Milky Way and Andromeda already share a merged gravity well explaining their future merger.

Φ(x,y)  =  i  GMidieκdi,di=(xxi)2+(yyi)2 \Phi(x,y)\;=\;-\sum_{i}\;\frac{G\,M_i}{d_i}\,\mathrm{e}^{\kappa\,d_i}, \qquad d_i=\sqrt{(x-x_i)^2+(y-y_i)^2}\,
grid (px)
half-span (kpc)
κ (per kpc)
(≈ 2.59e-1 m⁻¹)
Try κ ≈ 0.005–0.02/kpc
1000–1200 kpc shows full MW–M31 bridge.

Supercluster Flow (2D)

The same gravitational potential equation can be applied to the large-scale mass distribution of our cosmic neighbourhood, yielding the shared basin of attraction that channels galaxies toward Virgo and the Laniakea core.

Φ(x,y)  =  i  GMidieκdi,di=(xxi)2+(yyi)2 \Phi(x,y)\;=\;-\sum_{i}\;\frac{G\,M_i}{d_i}\,\mathrm{e}^{\kappa\,d_i}, \qquad d_i=\sqrt{(x-x_i)^2+(y-y_i)^2}\,
grid (px):
half-span (Mpc):
κ (per Mpc):
arrow step (px):

Φ uses the same κ factor as rotation/lensing: higher κ deepens wells over large d.

Flow arrows trace **infall** toward attractors (Laniakea core, Virgo, Shapley, etc.).

Use span ≈ **300–400 Mpc** to view the larger context; **120–200 Mpc** for Local Group + Virgo.

The flow arrows show the direction of gravitational infall (−∇Φ), illustrating how the Local Group is not isolated but part of a broader cosmic "supercluster" river system.

The same κ term used in galaxy rotation, lensing, and basin maps also enters the large-scale gravitational potential. When averaged over cosmological distances—dominated by voids rather than dense structures—it produces a small net positive contribution to the integrated potential: an emergent large-scale acceleration.

Φ(r)  =  GMreκr \Phi(r) \;=\; -\frac{GM}{r}\,e^{\kappa r}

For large radii, expanding the exponential gives an effective acceleration

a(r)  =  Φ    GMr2(1+κr), a(r) \;=\; -\nabla\Phi \;\approx\; -\frac{GM}{r^2}\,\bigl(1 + \kappa r\bigr),

so that κ contributes a small outward term proportional to κ on large scales. When this contribution is averaged over the cosmic web, it acts in the same direction as a cosmological constant, but arises from structure rather than vacuum energy.

Effective acceleration term in the Friedmann equation

In a homogeneous background, the large-scale effect of κ can be summarised as an additional acceleration term in the Friedmann equation:

a¨a=4πG3ρeff  +  Aκ, \frac{\ddot{a}}{a} = -\frac{4\pi G}{3}\,\rho_{\text{eff}} \;+\; \mathcal{A}_\kappa,

where 𝒜κ is an effective contribution generated by the large-scale κ field. For a representative background value κ₀ ≈ 2.6×10−26 m−1 (from supercluster flows), the associated acceleration scale 𝒜κ is of the same order of magnitude as the late-time acceleration usually attributed to Λ in ΛCDM.

In this view, the observed cosmic acceleration emerges from the cumulative effect of structure-dependent curvature, not from a fundamental vacuum energy term.

The Hubble Tension

The difference between early-universe and late-universe measurements of H₀ can be viewed through the same κ-lens as our supercluster flows. Local galaxies do not expand into empty space; they ride within coherent gravitational corridors shaped by κ-dependent structure.

Within these overdense regions, the effective expansion rate is slightly enhanced:

H0(κ)    H0(CMB)(1+βκrlocal) H_0^{(\kappa)} \;\simeq\; H_0^{(\text{CMB})} \left(1 + \beta\,\kappa\,r_{\text{local}}\right)

where β ≈ 1–2 parameterises how strongly local κ-dependent flows couple to the global expansion.

For a representative κ ≈ 8×10−4 Mpc−1(corresponding to κ₀ ≈ 2.6×10−26 m−1) and rlocal ≈ 100 Mpc with β ≈ 1.1:

H₀(κ) ≈ 67 × (1 + 0.09) ≈ 73 km s⁻¹ Mpc⁻¹

This illustrates that the same κ–driven structural acceleration that shapes basin and supercluster flows can naturally generate a 5–10% enhancement in the locally inferred H₀, comparable to the Planck–SH₀ES tension.

Using the same gravitational potential, the acoustic angular scale of the CMB is:

θ=rs(z)DA(z),πθ. \theta_\star=\frac{r_s(z_\star)}{D_A(z_\star)}, \quad \ell_\star \simeq \frac{\pi}{\theta_\star}.

θ_* ≈ 144.6 Mpc / 13.9 Gpc ≈ 0.0104 rad ≈ 0.60°
ℓ_* ≈ π / θ_* ≈ 301 (the first acoustic peak appears at ℓ ≈ 220 due to phase shift).

Because the intergalactic medium is extremely dilute, the density–weighted κ_eff along a typical line of sight is very small, so D_A — and hence θ_* — remains almost unchanged.

κeff=1L0Lk0 ⁣(ρ(s)ρ0)ads,DA(κ)DAexp ⁣(12κeffL). \kappa_{\text{eff}} = \frac{1}{L}\int_0^L k_0\!\left(\frac{\rho(s)}{\rho_0}\right)^{a}\,ds, \qquad D_A^{(\kappa)} \approx D_A\,\exp\!\Big(\tfrac12\,\kappa_{\text{eff}}L\Big).

With a void–dominated line of sight:
κ_eff ≈ 3×10⁻²⁹ m⁻¹ and L ≈ 4.3×10²⁶ m → ½ κ_eff L ≈ 0.0065, so D_A^(κ) / D_A ≈ exp(0.0065) ≈ **1.0065** (≈ +0.65%).

Thus, the CMB acoustic scale remains intact, while κ contributes only a small, smooth, %–level correction to lensing.

ακ(b)=αGR(b)eκb/2\alpha_\kappa(b)=\alpha_{\rm GR}(b)\,e^{\kappa b/2}

Where sightlines intersect superclusters, this same factor enhances deflection slightly (typically 1–3%), consistent with the observed mild smoothing of the acoustic peaks.

Post-Newtonian Limit: GR Locally, κ₀ as a Small Correction

On Solar-System scales, any modification of gravity must reduce to the standard post-Newtonian form tested by planetary orbits and light deflection. In the κ–r framework this can be achieved by treating the κ-response as a very small correction to the time–time component of the metric.

gtt;=;eκ(r)r,κ(r)κ02GMc2r2+,U=GMc2rg_{tt} ;=; -e^{\kappa(r)\,r},\quad \kappa(r) \approx \kappa_{0} - \dfrac{2GM}{c^{2}r^{2}} + \cdots,\quad U = \dfrac{GM}{c^{2}r}
ds2(12U+κ0r)c2dt2  +  (1+2U)(dr2+r2dΩ2)  +  O(c4)ds^{2} \simeq -\Big(1 - 2U + \kappa_{0} r\Big)c^{2}dt^{2} \; + \; \Big(1 + 2U\Big)\,(dr^{2}+r^{2}d\Omega^{2}) \; + \; \mathcal{O}(c^{-4})
γ1,β1,κ0rSolar1\gamma \simeq 1,\quad \beta \simeq 1,\quad \kappa_{0} r_{\text{Solar}} \ll 1
κ02.6×1026 m1 \kappa_{0} \approx 2.6 \times 10^{-26}\ \text{m}^{-1}

With κ₀ at this level, the extra κ₀ r term in gtt is completely negligible on Solar-System scales, so the standard post-Newtonian parameters remain indistinguishable from their GR values within current bounds, while still allowing a small cumulative effect from κ₀ on cosmological scales.

Mercury: the Famous 43″/Century Test

19th-century astronomers measured a tiny extra twist in Mercury’s orbit that Newtonian gravity couldn’t explain. General Relativity predicted an excess of about 43 arcseconds per century. The κ–r geometry matches the same result locally (with no extra parameters), and any cosmological bias from κ0 is far below detectability.

κ0=0 → pure local geometry (GR). Non-zero adds an (undetectably small) cosmological bias.
42.996″ / century
Δϕ=42.996arcsec/century\Delta\phi = 42.996\,\mathrm{arcsec/century}
QuantitySymbol / FormulaValue
Semi-major axisa = rp/(1−e)0.387073 AU   (5.791e+10 m)
Orbits per century36525 / 87.9691415.2
Per-orbit GR precession
ΔϕGR=6πGMc2a(1e2)\Delta\phi_{\rm GR} = \dfrac{6\pi GM}{c^{2} a (1-e^{2})}
0.10355″ / orbit
GR per centuryΔφGR × (orbits/century)42.996″ / century
κ0 correction× (1 + κ0 a)1.00000
Predicted precession42.996″ / century

Observed excess (over Newtonian/perturbative precession): ≈ 43.0″/century. With κ0=0 this panel reproduces the GR value. For κ0≈2.6×10⁻²⁶ m⁻¹, the additional shift is ~10⁻⁴″/century — below current detectability.

Interstellar Visitors in the Current κ(r) Model

The κ–response used for galaxies and clusters can also be applied to small bodies passing through the inner Solar System. In the present model,

κ(r)=κ0+kv(v/r1012s1)3(ρρ0)1/2, \kappa(r) = \kappa_{0} + k_{v} \left(\frac{\partial v / \partial r}{10^{-12}\,\mathrm{s}^{-1}}\right)^{3} \left(\frac{\rho}{\rho_{0}}\right)^{1/2},

with κ₀ ≈ 2.6×10⁻²⁶ m⁻¹, kᵥ ≈ 5×10⁻²⁶ m⁻¹ and ρ₀ = 1600 kg m⁻³. For a Keplerian profile v(r) = √(GM⊙/r) one has |∂v/∂r| = v/(2r), and for a representative interplanetary density ρ(1 AU) ≈ 10⁻¹⁹ kg m⁻³ scaling as 1/r², the resulting κ(r) at the heliocentric distances of 1I/ʻOumuamua, 2I/Borisov and 3I/Atlas is:

κr    {3×108(r0.25 AU, 1I)1.5×1011(r1.36 AU, 3I)2.6×1012(r2.0 AU, 2I) \kappa r \;\approx\; \begin{cases} 3\times 10^{-8} & (r \approx 0.25\ \text{AU, 1I})\\[4pt] 1.5\times 10^{-11} & (r \approx 1.36\ \text{AU, 3I})\\[4pt] 2.6\times 10^{-12} & (r \approx 2.0\ \text{AU, 2I}) \end{cases}

The effective solar acceleration is then

aeff(r)=GMr2eκ(r)r    aNewton(r)(1+κ(r)r), a_{\text{eff}}(r) = \frac{GM_\odot}{r^{2}}\, e^{\kappa(r)\,r} \;\simeq\; a_{\text{Newton}}(r)\,\bigl(1 + \kappa(r)\,r\bigr),

so that the fractional correction Δa/a is at most ≈ 3×10⁻⁸ for ʻOumuamua and much smaller for 2I/Borisov and 3I/Atlas. In the current κ(r) model the impact of κ on inner–Solar–System dynamics is thereforewell below present observational precision, while remaining large enough on galactic and cluster scales to account for rotation curves and lensing without dark matter.

Gravitational Waves in a κ–r Universe

Gravitational waves are one of our sharpest tests of gravity. In the κ–r geometry, present–day signals from neutron star and black hole mergers are indistinguishable from GR, while the same curvature response predicts enhanced primordial waves in the very early universe.

Φeff(r)=GMreκ(r)r\Phi_{\text{eff}}(r) = -\dfrac{GM}{r}\,e^{\kappa(r)\,r}
heff    hGReκ(r)rh_{\text{eff}} \;\propto\; h_{\text{GR}}\,e^{\kappa(r)\,r}
For κr1:eκr1+κr    heffhGR\text{For } \kappa r \ll 1:\quad e^{\kappa r} \simeq 1 + \kappa r \;\Rightarrow\; h_{\text{eff}} \simeq h_{\text{GR}}

Local mergers: GR recovered

Neutron–star and black–hole binaries live in regions where κ r ≪ 1, so the exponential factor is essentially unity.

Phase evolution, chirp mass and waveform shape reduce to standard GR:

gμν(κ)gμνGR(Solar System / stellar densities)g_{\mu\nu}^{(\kappa)} \simeq g_{\mu\nu}^{\rm GR} \quad (\text{Solar System / stellar densities})

For GW170817–like systems, the κ–r model reproduces a strain of h ∼ 4×10⁻²¹, matching LIGO/Virgo observations.

Early universe: enhanced primordial waves

In the very early universe, densities and velocity gradients drive κ(r) to much larger values, so κ r ≳ 1.

The same factor that is negligible today becomes important:

hprim    hGR,primeκearlyrh_{\text{prim}} \;\propto\; h_{\text{GR,prim}}\,e^{\kappa_{\text{early}} r}

This predicts a modest enhancement of the primordial gravitational–wave background and associated CMB B–modes, providing a clean target for future missions.

Today's detectors therefore see GR–exact waveforms, while the earliest gravitational waves are subtly reshaped by κ(r). The κ–r model passes current tests and makes falsifiable predictions for primordial signals.

Supermassive Black Holes: Born Heavy

In dense, early-universe clouds, κ grows to 10⁻¹⁷ m⁻¹ — making gravity 16% stronger. Collapse accelerates. Accretion explodes. A 10⁹ M⊙ black hole forms in under 10 million years.

κ5×1017 m1,eκr1.16,tcollapse0.93tff \kappa \sim 5 \times 10^{-17}\ \text{m}^{-1},\quad e^{\kappa r} \sim 1.16,\quad t_{\text{collapse}} \sim 0.93 \, t_{\text{ff}}

TOV Baseball: A Neutron Star in Your Hand

Imagine a fully loaded baseball diamond of neutron stars — four 1.4 M⊙ stars at the corners, 100,000 meters apart. Each packed with ρ ≈ 6.0 × 10¹⁷ kg/m³.

κ5×1017 m1,eκr1.16 \kappa \approx 5 \times 10^{-17}\ \text{m}^{-1},\quad e^{\kappa r} \approx 1.16

The central acceleration jumps from 0.85 m/s² to 0.99 m/s² — enough to trigger Schwarzschild collapse in under 1.5 km.

This shows how κ amplifies collapse in dense environments — the same mechanism that drives rapid SMBH formation in the early universe.

See PDF Section 3.4.1: "The TOV Baseball"

Descent: The Quantum Limit

If κ encodes structure at every scale, where does that structure end?
What happens when r → ℓ_P — the quantum domain where mass and weight separate?

Φ(r)=GMreκr,limrPΦ(r)=GMr \Phi(r) = -\frac{GM}{r}\,e^{\kappa r}, \quad \lim_{r \to \ell_P} \Phi(r) = -\frac{GM}{r}

At Planck scales, κ loses leverage. Curvature decouples from structure.
The exponential vanishes, restoring the unweighted Newtonian (and GR) potential.

Eκ=mc2eκrEκmc2(as r0) E_\kappa = m c^2\,e^{\kappa r} \quad \rightarrow \quad E_\kappa \to m c^2 \quad (\text{as } r \to 0)

Energy gain vanishes at small r — but seeds the first structure at larger scales.

Scale: r = 1.00e-35 m

κ r = 0.000000

eκr = 1.000000

Planck Scale

Φ(r)=GMr×1.0000 \Phi(r) = -\frac{GM}{r} \times 1.0000

κ → 0: Pure GR

The transition defines a natural cutoff: below it, mass is inertial; above it, it carries geometric weight.
See PDF §3.8: "Quantum Scale Indications"

Mass–Energy Equivalence in κ–Modified Gravity

The rest–energy relation E = mc² remains unchanged in the κ–model. Mass retains its inertial role. What changes is how energy couples to curvature. The effective gravitational mass acquires a scale–dependent weight through the factor exp(κ·r).

mgrav(r)=meκr m_{\text{grav}}(r) = m\,e^{\kappa r}

This introduces a distinction between inertial mass and gravitational mass without altering local special relativistic physics. At small radii, the exponential term approaches unity.

limr0mgrav=m,limr0Eκ=mc2 \lim_{r \to 0} m_{\text{grav}} = m,\qquad \lim_{r \to 0} E_{\kappa} = mc^{2}

At galactic and cluster scales, the κ-term enhances gravitational interactions by weighting energy according to local density and shear. At quantum scales, the weighting disappears, and the conventional mass–energy equivalence governs the dynamics.

See Appendix A.6: “Mass–Energy Equivalence Under κ(r)”

Appendix: Key Derivations

This appendix outlines the main steps behind the κ–modified gravity equations used in the text. Each derivation is shown in a compact, weak–field form suitable for galaxies, clusters, and large–scale structure.

1. Exponential potential from modified curvature

The starting point is an exponential f(R) action:

S=g[ReαR+16πGLm]d4x S = \int \sqrt{-g}\,\big[ R\,e^{\alpha R} + 16\pi G\,\mathcal{L}_m \big]\, d^4x

Varying this action with respect to the metric gμν gives the modified field equations:

f(R)Rμν12f(R)gμνμνf(R)+gμνf(R)=8πGTμν, f'(R) R_{\mu\nu} - \tfrac{1}{2} f(R) g_{\mu\nu} - \nabla_\mu \nabla_\nu f'(R) + g_{\mu\nu} \Box f'(R) = 8\pi G\,T_{\mu\nu},

where f(R) = R e\alpha R and

f(R)=ddR(ReαR)=eαR(1+αR). f'(R) = \frac{d}{dR}\big( R e^{\alpha R} \big) = e^{\alpha R}\,(1 + \alpha R).

In the weak–field regime relevant for galaxies and clusters, the curvature is small and |αR| ≪ 1. The exponential then admits the series expansion:

eαR1+αR. e^{\alpha R} \approx 1 + \alpha R.

To leading order, the corrections appear as small, R–dependent terms in the effective Poisson equation. Solving the modified field equations for a static, spherically symmetric mass M yields an effective potential that can be written in the form

Φeff(r)GMreκr, \Phi_{\rm eff}(r) \simeq -\,\frac{GM}{r}\,e^{\kappa r},

where κ collects the weak–field imprint of the exponential curvature term and depends on the local configuration of matter. The corresponding radial acceleration is

geff(r)=dΦeffdrGMr2eκr. g_{\rm eff}(r) = -\frac{d\Phi_{\rm eff}}{dr} \simeq \frac{GM}{r^2}\,e^{\kappa r}.

This is the universal κ–modified law used throughout the main text.

2. Orbital velocity and κ from rotation curves

For a test mass on a circular orbit of radius r around mass M, the centripetal acceleration is v² / r. Equating this to the κ–modified gravitational acceleration gives:

vκ2r=GMr2eκr. \frac{v_\kappa^2}{r} = \frac{GM}{r^2}\,e^{\kappa r}.

Solving for vκ:

vκ(r)=GMreκr/2. v_\kappa(r) = \sqrt{\frac{GM}{r}}\, e^{\kappa r / 2}.

The Newtonian prediction from baryonic mass alone is

vN(r)=GMr. v_N(r) = \sqrt{\frac{GM}{r}}.

The ratio between the observed orbital speed vobs(r) and the Newtonian prediction defines an empirical κ at radius r:

vobs(r)vN(r)=eκ(r)r/2. \frac{v_{\text{obs}}(r)}{v_N(r)} = e^{\kappa(r)\,r/2}.

Solving this relation for κ(r) gives:

κ(r)=2rln ⁣(vobs(r)vN(r)). \kappa(r) = \frac{2}{r} \ln\!\bigg( \frac{v_{\text{obs}}(r)}{v_N(r)} \bigg).

This expression is used to derive κ(r) directly from rotation curve data, without assuming any dark matter halo. The environmental model κ(r) in the main text is then fitted to these inferred κ values.

3. Gravitational lensing with κ

In standard General Relativity, the deflection angle for a light ray passing a mass M with impact parameter b is

αGR(b)=4GMc2b. \alpha_{\rm GR}(b) = \frac{4GM}{c^2 b}.

In the κ model, the same exponential correction that modifies the potential also modifies the lensing deflection. In the weak–field limit, the effective deflection angle can be written as

αeff(b)=αGR(b)eκb/2=(4GMc2b)eκb/2. \alpha_{\rm eff}(b) = \alpha_{\rm GR}(b)\, e^{\kappa b / 2} = \left( \frac{4GM}{c^2 b} \right) e^{\kappa b / 2}.

For κb ≪ 1, this reduces to

αeff(b)αGR(b)(1+12κb), \alpha_{\rm eff}(b) \approx \alpha_{\rm GR}(b)\, \big(1 + \tfrac{1}{2}\kappa b\big),

showing that κ introduces a small, scale–dependent enhancement to lensing without changing the underlying baryonic mass.

4. Environmental κ(r) from shear and density

The geometric origin suggests that κ should depend on local structure. A simple observationally–motivated form used in the main text is

κ(r)=κ0  +  kv(v/r1012s1)3(ρρ0)1/2. \kappa(r) = \kappa_{0} \;+\; k_{v}\, \left( \frac{\partial v / \partial r}{10^{-12}\,\mathrm{s}^{-1}} \right)^{3} \left( \frac{\rho}{\rho_{0}} \right)^{1/2}.

Here:

κ₀ — background curvature term
kᵥ — shear–response coefficient
∂v/∂r — local velocity gradient (shear)
ρ / ρ₀ — density relative to a fiducial scale

The cubic dependence on the velocity gradient emphasises regions with strong shear (for example, spiral arms or shocked gas in cluster mergers), while the square–root dependence on density captures the enhanced curvature in compressed structures relative to diffuse environments.

When κ(r) defined this way is inserted back into the expressions for vκ and αeff, the resulting predictions match observed rotation curves and lensing profiles across a wide range of systems using only baryonic matter.

5. Large–scale κ and an effective acceleration term

On very large scales, κ is dominated by the average properties of the cosmic web: voids, filaments, walls, and superclusters. In this regime, κ can be approximated by a slowly varying background value κ₀.

In a homogeneous background, the κ–modified gravitational response appears as an additive term in the acceleration equation,

a¨a=4πG3ρeff  +  A(κ0), \frac{\ddot{a}}{a} = -\frac{4\pi G}{3}\,\rho_{\rm eff} \;+\; \mathcal{A}(\kappa_0),

where 𝔄(κ₀) is an effective acceleration term built from the large–scale κ field. For suitable choices of κ₀ consistent with structure formation, this term can mimic a small, positive late–time acceleration similar in magnitude to the observed cosmological constant, without introducing a separate dark energy fluid.

The detailed identification of 𝔄(κ₀) with a specific Λ–like parameter depends on the averaging scheme and lies beyond the weak–field derivations used for galaxies and clusters, but the qualitative behaviour follows directly from the same κ–dependent correction to the potential.

6. Quantum limit of κ

We start from the κ–weighted potential used in the main text:

Phikappa(r)=fracGMr,ekappa(r),r \\Phi_{\\kappa}(r) = -\\frac{GM}{r}\\,e^{\\kappa(r)\\,r}

Here \\(\\kappa(r)\\) encodes the response of gravity to large–scale structure (background curvature, shear, and density). To understand the behaviour near the quantum limit, we examine \\(r \\to \\ell_P\\), where the Planck length \\(\\ell_P\\) is the characteristic scale below which classical structure is no longer resolved.

6.1 Small–r expansion of the exponential

For any finite \\(\\kappa(r)\\), the exponential admits a Taylor expansion around \\(r = 0\\):

ekappa(r),r=1+kappa(r),r+tfrac12,kappa(r)2r2+mathcalO(r3). e^{\\kappa(r)\\,r} = 1 + \\kappa(r)\\,r + \\tfrac{1}{2}\\,\\kappa(r)^{2} r^{2} + \\mathcal{O}(r^{3}).

Substituting into \(\Phi_{\kappa}(r)\) gives:

Phikappa(r)=fracGMrBigl[1+kappa(r),r+tfrac12,kappa(r)2r2+mathcalO(r3)Bigr] \\Phi_{\\kappa}(r) = -\\frac{GM}{r} \\Bigl[ 1 + \\kappa(r)\\,r + \\tfrac{1}{2}\\,\\kappa(r)^{2} r^{2} + \\mathcal{O}(r^{3}) \\Bigr]

Expanding term by term:

Phikappa(r)=fracGMr;  GM,kappa(r);  tfrac12,GM,kappa(r)2r;+  mathcalO(r2) \\Phi_{\\kappa}(r) = -\\frac{GM}{r} \\;-\; GM\\,\\kappa(r) \\;-\; \\tfrac{1}{2}\\,GM\\,\\kappa(r)^{2} r \\;+\; \\mathcal{O}(r^{2})

The leading term is the usual Newtonian potential \\(-GM/r\\). The \\(\\kappa\\)-dependent terms are finite or vanish as \\(r \\to 0\\), so the short–distance \\(1/r\\) structure of gravity is unchanged.

6.2 κ sourced by macroscopic structure

In the κ–r model, \\(\\kappa(r)\\) is an effective parameter built from coarse–grained structure:

kappa(r)=kappa0;+  kvleft(fracpartialv/partialr1012,mathrms1right)3left(fracrhorho0right)1/2 \\kappa(r) = \\kappa_{0} \\;+\; k_{v} \\left( \\frac{\\partial v / \\partial r}{10^{-12}\\,\\mathrm{s}^{-1}} \\right)^{3} \\left( \\frac{\\rho}{\\rho_{0}} \\right)^{1/2}

At Planck scales, matter distribution is effectively homogeneous and gradients vanish. Therefore:

limrtoellPkappa(r)=0 \\lim_{r \\to \\ell_P} \\kappa(r) = 0

And the κ–weighted potential reduces to:

limrtoellPPhikappa(r)=fracGMr \\lim_{r \\to \\ell_P} \\Phi_{\\kappa}(r) = -\\frac{GM}{r}

6.3 κ–weighted mass–energy

Ekappa(r)=mc2,ekappa(r)r E_{\\kappa}(r) = mc^{2}\\,e^{\\kappa(r) r}

Expanding for small r:

Ekappa(r)=mc2bigl[1+kappa(r)r+tfrac12kappa(r)2r2+mathcalO(r3)bigr] E_{\\kappa}(r) = mc^{2} \\bigl[ 1 + \\kappa(r) r + \\tfrac{1}{2} \\kappa(r)^{2} r^{2} + \\mathcal{O}(r^{3}) \\bigr]

Giving the limit:

limrtoellPEkappa(r)=mc2 \\lim_{r \\to \\ell_P} E_{\\kappa}(r) = mc^{2}

6.4 Interpretation

κ acts as a structural modifier: it vanishes when structure cannot be resolved (Planck scale) and grows when gradients, density contrasts, and shear appear on macroscopic scales.

Below the Planck scale, gravity reverts to its standard form. Above it, κ encodes geometric weight.

A.7 — Mass–Energy Equivalence Under κ(r)

In the κ–modified weak–field limit, the effective gravitational potential takes the form

Φκ(r)=GMreκr. \Phi_{\kappa}(r) = -\,\frac{GM}{r}\,e^{\kappa r}.

Differentiating gives the radial acceleration:

gκ(r)=GMr2eκr. g_{\kappa}(r) = \frac{GM}{r^{2}}\,e^{\kappa r}.

This may be interpreted as the usual Newtonian term multiplied by a scale–dependent gravitational weight. Writingmgrav = m · e^κ r reproduces the same force law.

F=mgravGMr2withmgrav(r)=meκr. F = m_{\text{grav}}\,\frac{GM}{r^{2}} \quad\text{with}\quad m_{\text{grav}}(r) = m\,e^{\kappa r}.

Inertial mass remains unchanged, so the rest–energy relationE = mc² holds exactly. The gravitational contribution to the energy, however, acquires the same weight:

Eκ(r)=mc2eκr. E_{\kappa}(r) = mc^{2}\,e^{\kappa r}.

At small radii, the weighting disappears and the standard expression is recovered:

limr0Eκ(r)=mc2. \lim_{r \to 0} E_{\kappa}(r) = mc^{2}.

This establishes a scale–dependent distinction between inertial and gravitational mass without altering local special–relativistic physics. Energy retains its inertial identity, while its gravitational influence varies with structure through κ(r).

Full derivations, MCMC fits, and code at: